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Abstract 


A mathematical model is developed and used to investigate the effects of axially distributed current collection on the performance of an 
anode-supported tubular solid oxide fuel cell (SOFC). The physical model considers coupled fluid and thermal transport within the tube, reactive 
porous-media flow within the electrodes, electrochemical charge-transfer, thermal transport within the MEA structure, and an electric circuit to 
represent the current-collection system. Four examples are used to illustrate the model and evaluate current-collection strategies. All the examples 
consider a particular Ni-YSZ | YSZ | LSM—YSZ tube, but with different current collection. The inlet fuel stream is a syngas mixture that results 
from steam reforming of dodecane. Results show that current collection can strongly affect local performance (i.e., composition, temperature, and 
current density profiles) along the length of the tube as well as overall performance (i.e., efficiency and utilization). 


© 2007 Elsevier B.V. All rights reserved. 
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1. Introduction 


Tubular cells and stacks have some advantages compared to 
planar cells and stacks. For example, ceramic-to-metal seals are 
usually confined to small areas at the ends of the tubes. Fuel- 
and air-distribution manifolds may be more straightforward in 
tubular stacks. Stresses caused by thermal cycling can be less 
troublesome in tubular systems because the membrane-electrode 
assembly (MEA) is less rigidly confined. However, current col- 
lection is much more challenging in a tubular cell, especially for 
anode current collection on the inside of the tube. This paper 
develops a physically based computational model that can be 
used to evaluate alternative current-collection strategies. 

Fig. 1 illustrates the nominal geometry for a tubular cell. This 
paper considers a porous cermet anode-supported tube, with the 
thin dense electrolyte and porous cathode applied to the outside 
of the tube. A typical tube wall may be fabricated as Ni-YSZ 
that is on the order of a millimeter thick. The dense electrolyte 
membrane is typically on the order of 10-20 um of YSZ. A typ- 
ical cathode may be 50 um of porous LSM-YSZ. As illustrated 
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in Fig. 1, the current collection is achieved by metallic wires 
that are bonded to the electrode surfaces. Additionally, a low- 
electrical-resistance bus is assumed to run axially the length of 
the tube. 

In principle, because the electrodes have electrical conduc- 
tivity, current could be collected axially through the electrode 
structures without any explicit current-collection structure. 
However, using typical electrode materials, the voltage drop 
would be far too high to be practical. The cathode especially is 
very thin and uses ceramic mixed ionic—electronic conducting 
(MIEC) materials that are relatively poor electronic conduc- 
tors. The Ni in an Ni-YSZ has high electronic conductivity 
and the structure is relatively thick. Nevertheless, a dedicated 
low-resistance current-collection significantly improves perfor- 
mance. 

The model in this paper considers fluid flow and heat- 
transfer within the tube, transport and catalytic chemistry 
within the porous electrodes, electrochemical charge-transfer, 
and heat-transfer within the MEA structure. An electrical net- 
work describes the current-collection system. Circumferential 
current-collection wires are attached to the electrodes at speci- 
fied positions along the length of the tube. Usually, because of 
higher electrical-resistance in thin MIEC cathode materials, the 
cathode attachment is more frequent than the anode attachment. 
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Fig. 1. Sketch of an anode-supported tubular SOFC. 


An electrical-resistance is assigned to each attachment point, 
which may represent interface or bonding resistance. The 
attachment resistance may also account for the resistance to cir- 
cumferential current flow in the wires leading to an axial bus. 
Each segment of the axial bus, between attachment points of the 
circumferential wires, is also assigned a resistance. The elec- 
trodes also have electrical-resistance. The model accounts for 
the local electric-potential profiles within the electrodes. 

There is recent literature concerning the fabrication and per- 
formance of tubular cells [1-4]. These works are concerned 
primarily with ceramics materials processing. A number of 
recent modeling studies have also been developed to investigate 
the performance of tubular cell and systems [5-15]. 

Nagata et al. [5] and Aguiar et al. [6,7] solve energy and mass 
balances to represent one-dimensional steady-state reacting- 
flow for single-cell tubular SOFCs. Other groups have developed 
two-dimensional axisymmetric models [8-15]. These models 
predict pressure, temperature and species concentration dis- 
tributions. In some cases circumferential current flow results 
in loss of full axisymmetry. Nevertheless, chemistry and heat 
release are based on circumferentially averaged electric poten- 
tials. 

Charge-transfer chemistry is usually represented in Butler— 
Volmer form. Li et al. [8-11] assume that H2 is the only 
electrochemically active species. They also assume that CH4 
reforming and water-gas-shift reactions are in equilibrium. 
Suwanwarangkul et al. [14] assume that Hz and CO can both be 
electrochemically oxidized at the anode-electrolyte interface. 
They use a global rate expression for the water-gas-shift reac- 
tion. Nishino et al. [15] also considered that both Hz and CO are 
electrochemically active. They represent CH4 steam reforming 
and water-gas-shift reactions by the global kinetics. Campanari 
and Iora [12] assume that hydrogen is the only electrochemi- 
cal active species. They use a global kinetic expression for CH4 
steam reforming and assume that the water-gas-shift reaction 
remains equilibrated. 


Boersma et al. [16] develop an electric-network model for 
a single-tube fuel cells to describe the axial electric-potential 
distribution. The model incorporates ohmic losses across the 
MEA and along the length of the tube. For the particular system 
studied, the model predicts that the active tube length should 
be no longer than 10 cm due to the ohmic losses along the tube 
length, and the electrolyte should be less than 20 wm thick to 
achieve high power density. Suzuki et al. [17,18] have compared 
the current-collection efficiency for micro-tubular SOFCs, using 
either single-terminal or double-terminal anode current collec- 
tion. They report that double-terminal anode current collector is 
much more efficient than single-terminal current collector. 

The present paper considers the coupled effects of the current 
collection on fluid flow and chemistry, electric-potential distri- 
bution, and thermal behavior of the MEA structure and gas flow. 
The current-collection system, which directly affects the local 
cell potential, has an important impact on the overall cell per- 
formance. Thus, the current-collection system is an important 
aspect of tubular cell design and predictive models can greatly 
assist this process. 


2. Mathematical models 


The mathematical model in this paper builds on the model 
developed by Zhu et al. for planar cells [19]. In addition to 
considering a tubular geometry, the present model predicts the 
thermal response of the flow and the tube wall. The present 
model also incorporates the electrical behavior within the MEA 
and the current-collection system. The internal reforming chem- 
istry and the charge-transfer models are the same as those in Zhu 
et al. [19]. 


2.1. Flow within the tube 


The gas-phase flow in the tube is treated as a one-dimensional 
laminar flow, neglecting the spatial variations in the radial direc- 
tion [20]. For syngas or methane fuels, homogeneous chemical 
kinetics are negligible for temperatures below around 900° C 
[21]. Using a plug-flow approximation, the mass and momentum 
conservation equations are summarized as 
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These equations are written in transient form, with the time ¢ and 
the axial coordinate x being the independent variables. Depen- 
dent variables include the mass density p, the mean velocity 
u, the pressure p, and the species mass fractions Yg. Geomet- 
ric parameters include the hydrodynamic perimeter Ph and the 
cross-sectional flow area Ac. Assuming an inner tube diameter 
D, Py = xD and A, = aD? /4. 
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The density is determined from the pressure and composition 
using an ideal-gas equation of state 


pw Pp 
= = 4 
P= RT Ka ; (4) 
RTX Yx/ We 
k=1 


where R is the universal gas constant, T the temperature, W 
the mean molecular weight, and W% are the species molecular 
weights. Kg is number of gas-phase species. 

As a result of internal reforming and charge-transfer chem- 
istry, there are species fluxes between the porous-anode structure 
and the flow within the tube. These molar fluxes are represented 
as a As discussed in the following section, these fluxes are 
determined from the solution of the radial reactive porous-media 
problem within the MEA. 

In addition to axial convective transport, there may be axial 
diffusion of species along the length of the tube. The axial dif- 
fusive mass fluxes are represented as 

We OX, 
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The mixture-average diffusion coefficients are calculated as, 
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where Xz are the mole fractions, and Dx are the binary diffusion 
coefficients [20]. 

The wall shear stress ty, is represented in terms of a friction 
factor f = 2t,,/pu*. For a circular channel, the friction factor 
can be calculated based on the laminar, fully-developed flows as 
Re f = 16. The local Reynolds number is defined in terms of 
the mean velocity u and the hydraulic diameter as Re = puD/u 
where u is the mixture viscosity. 


2.2. Porous-electrode model 


Because the electrodes are relatively thin and the axial gra- 
dients are low, axial and azimuthal species transport may be 
neglected within the porous anode structure. Radial transport 
and chemistry can be represented in one-dimensional form with 
overall and species continuity equations as 


APpYn) aM 


k = Aoi, We, 7 

at + or sSk Wk (7) 
K K 

agp) AGO A. 

eit’ X =) Acs; We. 8 

at t2 ar - sSk Wk (8) 


In these equations time f and radius r are the independent vari- 
ables. The electrode porosity is @ and Ag is the specific surface 
area of the active catalysts (i.e., active surface area per unit vol- 
ume of the porous media). The molar species production rates 
by heterogeneous reaction are represented as 5,. 


The gas-phase species mass fluxes hg through the pore struc- 
ture are determined from the Dusty-Gas Model (DGM) [22,19], 
which can be written as an implicit relationship among the 
gas-phase species molar fluxes a molar concentration [Xx], 
concentrations gradients, and the pressure gradient as, 
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where [Xr] = p/RT is the total molar concentration, B the per- 
meability, and jz is the mixture viscosity. The radial mass fluxes 
are related simply to the molar fluxes as Phe = WJ a The effec- 
tive binary and Knudsen diffusion coefficients Df, and Di xn are 
written as 
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Die = Ê Du, Dkn 
The binary diffusion coefficients Dge and the mixture viscosities 
u are determined from kinetic theory [20]. Knudsen diffusion 
represents mass transport assisted by gas—solid collisions. The 
Knudsen diffusion coefficients depend on physical character- 
istics of the porous media, including porosity ¢, average pore 
radius rp, primary particle diameter dp, and tortuosity t. The 
permeability can be determined from the Kozeny—Carman rela- 
tionship as 


3 72 
bd, 
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Further details of the DGM and its numerical implementation 
can be found in Zhu et al. [19]. 

Heterogeneous chemical kinetics proceeds on the catalyst 
surfaces within the porous anode. Molar production rates by 
the heterogeneous reactions 5, are functions of temperature, 
gas concentrations, and surface-species coverages. The reaction 
mechanism used in the present model incorporates steam and 
dry reforming as well as partial oxidation [23]. This mecha- 
nism considers 42 reactions among 6 gas-phase species and 12 
surface-adsorbed species. Temporal variations of the site cov- 
erages 6; of the K, surface-adsorbed species are represented 
as 
Stal peleni a2) 

dt r 
where T is the available site density. At steady-state, 5, = 0. 

Boundary conditions are needed to solve Eqs. (7) and (8) 
within the porous electrodes. At the interfaces between the 
porous electrode and the external gas flow (i.e., fuel within 
the tube or air outside the tube) the gas-phase composition is 
assumed to be that within the either the fuel or air flow. At 
the electrode-electrolyte interfaces, species fluxes are estab- 
lished by the electrochemical charge-transfer reactions within 
thin electrode-electrolyte three-phase boundary (TPB) regions 
as 

Wkyk 
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where ie is the local current density, F the Faraday’s constant, 
vz the stoichiometric coefficients in a global charge-transfer 
reaction and nę is the number of electrons transferred. For 
hydrogen oxidation at the anode-electrolyte interface (i.e., 
Hp + O% @ H20 + 2e7), vp, = —1, vmo = 1, and ne = 2. 
For oxygen reduction at the cathode-electrolyte interface (i.e., 
O2 + 4e— = 207-), uo, = 1 and ne = 4. 

The gas-flow equations and the porous-media equations are 
intimately coupled. The species molar fluxes between the surface 
of the porous electrode and the adjoining gas flow (i.e., the JM 
terms in Eqs. (1) and (3)) are evaluated from the porous-media 
solution at the interface. The gas-flow composition (determined 
from Eqs. (1) and (3)) serves as boundary conditions for the 
reactive porous-media problem (i.e., Eqs. (7) and (8)). 


2.3. Electrochemistry 


The operating cell potential between anode and cathode at the 
dense-electrolyte interfaces Ece may be represented in terms of 
overpotentials as 


Ece = Erev — lNact,c(ée)I, (14) 


where Erey is the reversible potential across the dense elec- 
trolyte membrane (i.e., not the reversible potential based on 
the electrochemical potentials of the gases in the anode and 
cathode compartments). The reversible potential E;ey, the cell 
potential Ece, and the current density i, all vary continuously 
along the length of the tube. The activation overpotentials at the 
anode-electrolyte and cathode-electrolyte interfaces are writ- 
ten as act and Nactc, respectively. The ohmic overpotential 
is written as Nohm. Because gas-phase transport and chemistry 
within the electrode structure are modeled fully in this formu- 
lation, there is no explicit need to incorporate a “concentration 
overpotential.” 

The current density is written in terms of the Butler—-Volmer 
equation as 


TOE Oa FMact Ae Fact 
ie = ig |exp RT exp RT ; (15) 


where ig is the exchange current density, F is the Faraday con- 
stant, and a, and a, are the anodic and cathodic symmetry 


Nact,a(le) — Nohm(ée) — 


parameters. 
The exchange current density at the anode is given as [19] 
( )" omo) 
b=, P/P) (Pmo (16) 


1+ (pup / Pi.) 


The parameter PH depends on hydrogen adsorption/desorption 
rates. The dependence of the exchange current density ig on the 
temperature can be expressed as, 


Ey, f1 1 
atari an 


where Ey, is the apparent activation energy for the electro- 
chemical oxidation of H2, and the parameter i*, m IS assigned 
empirically to fit measured polarization data at a reference tem- 
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perature Tef. Similarly, at the cathode, the exchange current 
density is expressed as 


J14 
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where 
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Eo, is the apparent activation energy for the electrochem- 
ical reduction of O2. The parameter Po, depends on O2 
adsorption/desorption rates and the parameter i* efO, 1S assigned 
empirically. All the partial pressures in Eqs. (16) and (18) are 
measured in atmospheres. 


2.4. Tube flow energy balance 


This model considers three temperatures: the gas tempera- 
ture within the anode tube T}, the tube-wall temperature Tm, and 
the external air temperature Te. The thermal balance among the 
fuel-cell components includes convective heat-transfer between 
flow in the tube and the MEA structure, convective heat-transfer 
between the outside air and MEA structure, heat transport asso- 
ciated with species transport between the electrodes and the 
adjoining gas, and axial heat conduction within the tube wall. 
Source terms appear in the energy balance for the tube wall 
(i.e., the MEA) to describe heat release associated with thermal 
reactions within the porous electrodes, electrochemical charge- 
transfer reactions, and ohmic heat generation. 

Energy conservation for gas flow within the tube can be 
expressed as 


ðE a P 
+ gh +), 


20 
ot ox G0) 


where E = pe and e is the specific internal energy of the gas. 
The heat flux g, which represents both axial heat conduction and 
gas-phase species transport, is written as 


Kg 


oT ; 
q = —A— + X (PYku + jih, 21) 
ox | 


where hg are the species specific enthalpies and À is the mixture 
thermal conductivity. 

The convective heat flux between the tube flow and the MEA 
structure is represented as 


q! = heony(T — Tm), (22) 


where /icony is a convective heat-transfer coefficient. Within 
the tube, Acony may be evaluated using a conventional Nusselt- 
number correlation. 

The second term on the right-hand side of Eq. (20) (i.e., gy 
represents the energy transport associated with mass transfer 
between the channel flow and the porous anode structure (i.e., 
tube wall). Because all the heat release (i.e., resulting from ther- 
mal and electrochemical reactions and ohmic heat resistance) is 
assigned to the MEA structure, qM must be subtracted from the 
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conservative form of the tube-flow energy equation (Eq. (20)). 
In the conservative form, the tube-flow energy equation implic- 
itly contains the heat released in the porous electrode through 
the species continuity equations. 

The heat flux associated with species transport between the 
MEA structure and the tube flow qM can be written as 


= $ I Webe(Ta) + X Wehr Tin) (23) 


JM>0 JM<0 


If species transfer from the tube into the MEA structure, 
enthalpies of the gas-phase species hg are evaluated at the tube 
gas temperature T,, otherwise, hę are evaluated at the MEA tem- 
perature Tm. Note that af > 0 indicates radially outward flux. 
Heat flux between the outside of the tube (cathode) and the air 
can be written as 


= — So IM Webe(Te) — So I Wihk(Ta). (24) 


IM <0 IM>0 


In this expression, the is the molar flux at the exterior surface 
of the cathode. In this case, pa > 0 means gas is leaving the 
cathode and entering the surrounding air. 

The heat-transfer coefficient within the tube Agony can be 
evaluated based on a Nusselt-number correlation as, 

heony D 
Nu = conv ; (25) 
À 

For the laminar flow in a circular channel with constant wall 
temperature, the Nusselt number can be calculated as [24,25], 


1000 \ 70488 57.2 
Nu = 3.657 + ost) exp (-22) : (26) 


where the Graetz number is defined as Gz = Re Pr D/x. 


2.5. MEA energy balance 


The MEA is assumed to have a uniform temperature radi- 
ally, but varying axially. Further, the solid materials and the gas 
within the pore spaces are at a common temperature Tm. Energy 
conservation for the MEA structure may be written as 


Em dqm ; Ph T Po T 
——_+— = + — qa + g 27 
at ag maT ee Am" oO 


where Enp is the total energy per unit volume of the MEA struc- 
ture, P, the outside perimeter of the tube and Am = nro — rž) 
is the cross-sectional area of the tube wall. The net volumet- 
ric heat source within the MEA structure, due to the thermal 
and electrochemical reactions and the ohmic resistance, may be 
written as 

mea = a q™ + Po X 


Pe E (28) 
i : 

An de An e£cell 

In this equation Pe is the perimeter of the dense electrolyte. Axial 

conduction within the MEA structure is written as 


Tm 
qm = im’ (29) 


where Àm is the effective thermal conductivity of the composite 
MEA structure. 

Eq. (27) is second order and requires boundary conditions at 
both ends of the tube. The results discussed later in this paper 
assume that the tube wall is insulated at both ends. 


3. Electrical network 


Electrons in the cathode participate in charge-transfer reac- 
tions to reduce the oxygen, and the electrons released by 
electrochemical fuel oxidation are delivered into the anode. The 
role of the current-collection architecture is to distribute elec- 
trons with minimal ohmic losses. Referring again to Fig. 1, we 
assume that the terminal voltage (i.e., cell potential) is fixed at 
one end of the tube. The current-collection system is formed by 
circumferentially wrapped high-conductivity wires that contact 
the electrodes and axial wires (buses) that carry current to the ter- 
minals. The electrode materials are also electrically conductive. 
However, it is presumed that the electrical conductivity within 
the electrodes (especially the cathode) is lower compared to the 
current-collection system. 

The net electrical-resistance (i.e., electrodes and current col- 
lection) serves to increase the local operating potential Ecen(x) 
along the length of the tube. This, in turn, limits the cur- 
rent that can be generated electrochemically. If the resistance 
is sufficiently great as to cause the local cell potential to 
reach the local open-circuit potential at some position x (i.e., 
Ecen(x) = Erey(x)), the remainder of the cell cannot produce 
electric current. In addition to primarily electrical aspects, there 
are potentially important consequences on thermal behavior 
and overall cell efficiency. Thus, careful design of the current- 
collection system to reduce net electrical-resistance is a very 
important aspect of tubular SOFC design. 

Fig. 2 illustrates an electric-circuit model that represents a 
single tubular SOFC cell. The cathode bus is represented as a 
series of resistances Re. This bus is connected to the cathode 


itself via an array of attachments that have resistances ee There 
is also axial electrical conduction within the cathode structure, 
which is represented through its effective electrical conductivity 
o£. In this circuit, the dense electrolyte behaves as an axially dis- 
tributed current source (much like a locally distributed battery). 
The anode side of the network is constructed in a way that is 
analogous to the cathode. Electrons flow from the porous anode 
structure through discrete attachment points into the anode bus, 
and ultimately to the cell terminal. As with the cathode, the anode 
structure has electrical conductivity of, and the attachment and 
bus have resistances RA and RB, , respectively. Because the 
anode and the cathode materials may have very different electric 
conductivities, the number of attachments and the resistances of 
the interconnect bus may be very different for the anode and the 
cathode current collection. 

The anode side of the electrical network may be represented 
with charge-conservation equations as 


ð e Ba Palxi) — DR, 
= (o: #2) E m EH), GO) 
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Anode 


Electrolyte 
Cathode 


i+l 


i+l 


= 


In these equations, ®, is the anode electric potential and DR, 
is the electric potential at the attachment node i of the anode 
interconnect bus. The Dirac delta function is ô, and n‘ is the 
number of attachments between the anode and the anode inter- 
connect bus. Eq. (30) is a balance equation for the local electric 
potential a(x) within the porous anode structure. The right- 
hand side has local source terms for the current production 
ie(x) and the current provided from the bus through the ns dis- 
crete attachment points. Eq. (31) represents the charge-balance 
at each attachment point i. Without capacitance in the circuit, 
the net rate of charge accumulation at each node must van- 
ish. This means that the sum of the currents at a node must 
vanish. 


Charge balance on the cathode side of the network is repre- 
sented as 


®. i — B, 
1 (o: =) aant S eigen), e 


ox Ox 4 RS 
B B B B B 
Poi = Poi- Poi = Ooi Poi Z Pelxi) =0 (33) 
RB RE tik T 
Gi c,i+1 ci 


In these equations, ®, is the cathode electric potential and pB, 
is the electric potential at the attachment node i of the cathode 
interconnect bus. 

The local electric-potential difference across the MEA struc- 
ture can be calculated as the potential difference between the 
cathode and the anode as 


Ecen lx) = P(x) — alx). (34) 


Assuming that the output operating cell voltage Fout is given, 

the terminal potential at the anode interconnect bus can be set 

to be zero as PBa = 0, and the terminal potential at the cathode 

interconnect bus can be set as the output operating cell potential 
Bi. 

as Poo = Egue 


4. Numerical algorithm 


The numerical solution of the coupled model is obtained by 
solving sequentially and iteratively three submodels: the fluid 
flow and heat transport along the channel length, the cell electric- 
potential distribution along the axial bus and the MEA structure, 
and the chemically and electrochemically reacting flows within 
the porous electrodes and electrolyte. At each iteration step, all 
the flow stream properties (i.e., gas composition Yg, tempera- 
tures Ta, Te, and Tm, and pressure p) and cell voltage Eve are 
first frozen. A one-dimensional MEA model is solved to obtain 
the electric current i, across the MEA structure, the heat fluxes 
gt and ge and the gas-phase species mass fluxes i between 
the flow streams and the electrodes. Then, by freezing the a 
ee and ie computed from the MEA model, the fluid flows and 
the temperatures along the channel are recalculated. The cell- 
voltage distribution along the channel is solved by freezing the 
current density ie through the MEA structure. 

Spatial derivatives in all the governing equations are dis- 
cretized using the finite-volume method. The transient problem 
of the MEA submodel is solved using a method-of-lines algo- 
rithm, with the time marching accomplished with the LIMEX 
software that is designed to solve systems of differential- 
algebraic Eq. [26]. The steady-state problem for the channel 
flow and the temperature is solved as a boundary-value problem 
using a hybrid-Newton method [27,20]. 


5. Results and discussions 


A particular anode-supported tube is used to illustrate the 
effects of the current-collection network on cell performance. 
The tube’s inner diameter is D = 0.8 cm and it is L = 25cm 
long. The anode is 900 um of porous Ni-YSZ, the electrolyte is 
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Table 1 

Physical and chemical parameters for the MEA structure 

Parameters Value 

Anode 
Thickness, La (wm) 900 
Porosity, @ 0.35 
Tortuosity, T 4.80 
Pore radius, Tp (um) 0.20 
Particle diameter, dp (um) 1.00 
Specific catalyst area, As (cm7!) 1080 
Apparent activation energy, Ey, (kJ mol`!) 120.0 
Reference temperature, Tef (K) 1073.0 
Exchange current factor, Def (A cm~?) 8.5 
Anodic symmetry factor, œa 1.5 
Cathodic symmetry factor, œe 0.5 

Cathode 
Thickness, Le (wm) 50 
Porosity, @ 0.35 
Tortuosity, t 4.00 
Pore radius, rp (wm) 0.25 
Particle diameter, dp (um) 1.25 
Apparent activation energy, Eo, (kJ mol`!) 130.0 
Reference temperature, Tef (K) 1073.0 
Exchange current factor, i760 (A cm~?) 2.4 
Anodic symmetry factor, œa f 1.5 
Cathodic symmetry factor, œe 0.5 


Electrolyte: ce, = ooT7! exp(— Ee / RT) 
Thickness, Le (wm) 20 
Activation energy, Ee, (J mol`!) 8.0E4 
Ion conductive pre-factor, oo (S cm7!) 3.6E5 


20 um YSZ, and the cathode is 50 um of porous LSM-YSZ. 
Table 1 lists other parameters that describe the MEA struc- 
ture. The apparent activation energies for the exchange current 
densities are taken from Nagata et al. [5]. The exchange-current- 
density factors iref, Hy and itp, are estimated by fitting the 
experimental data at the reference temperature Tref = 1073 K. 

The thermal conductivities of the anode, cathoed and 
electrolyte are taken to be 11.0, 6.23, and 2.7Wm!K7!, 
respectively [15]. The effective thermal conductivity of the MEA 
structure is estimated to be Am = 10.5 W m~! K7!. 

The effective anode and cathode electric conductivities 
(Scm7!) can be evaluated as functions of temperature in the 
form [28] 


e e 1 E° 
o = A'T exp ( =) (35) 
Values for the pre-factors A® and activation energies E° are listed 
in Table 2. Illustrative values of attachment and bus resistances 
also are shown in Table 2. The bus resistance R’ is stated per 
unit length. Thus, the net bus resistance between attachment 
points is RB = R’L. where L is the distance between attach- 
ments points. Assuming the buses are made of silver (electrical 
resistivity at 800°C, pe = 6.5 x 1076 Q cm), the bus resistances 
in Table 2 correspond to bus diameters of 3.2 mm for the cathode 
and 4.5 mm for the anode. 

In all cases the inlet fuel composition is 53.3% H2, 24.6% 
H20, 12.3% CO, 9.2% CO2, and 0.6% CH4. This is the equilib- 
rium output of a steam reformer operating on dodecane (C12H26) 


at 700°C, 1 atm, and a steam-carbon ratio of 2.5. The tube inlet 
velocity is 60cms~! and the inlet temperature is 750°C. The 
terminal voltage is fixed at Eout = 0.65 V. The model is devel- 
oped to consider the air depletion by solving the air flow through 
the air channel. But the air compartment for all the examples 
presented here is assumed to be flooded by the air such that 
the variations of species composition and temperature along 
the channel are negligible. The external air temperature and 
pressure are fixed to be 750° C and | atm. The heat-transfer 
coefficient between the cathode and the external air is taken to 
be 100 W m~? K7!. 

To explore the influence of current-collection strategies and 
to assist interpreting model results, four alternatives for current- 
collection systems are investigated. They are: 


Base case: The base case uses 50 cathode-attachment points 
and 10 anode-attachment points. 

Reduced cathode attachment: This case is the same as the base 
case, but with only 25 cathode attachment points. 

Without anode collection: This case assumes that the cathode 
electric potential is uniform (i.e. perfect cathode current collec- 
tion), and that all anode current collection is accomplished by 
axial conduction through the anode itself (i.e., no attachment 
points or anode bus). 

Uniform cell potential: This case assumes perfect current col- 
lection, such that the the cell potential Ece = Eout = 0.65 V 
everywhere along the tube. 


Fig. 3 shows a composite of computed solutions for the base 
case. The local cell potential and current density show a saw- 
tooth-like behavior, with the period set primarily by the cathode 
attachment points. The current (and power) density is nearly 
exactly out of phase with the cell potential. At the attachment 
points, where the cell potential is locally minimum, the current 
density is high. A low-frequency oscillation can also be seen, 
especially in the cell potential. This period is set by the less- 
frequent anode attachment. 

The bottom panel of Fig. 3 shows temperature and velocity 
profiles. Because of heat release due to chemistry, electrochem- 
istry, ohmic heat generation within the MEA, the tube-wall 
temperature increases in the entry regions where the current 
density is high. Relatively high thermal conductivity in the wall 
allows significant axial heat conduction. Fuel enters the tube at 


Table 2 

Electrical parameters for the MEA and current collection 

Parameters Value 

Anode 
Conductivity pre-factor, A§ (S K cm7!) 9.5 x 10° 
Conductivity activation energy, E$ (J mol`!) 9561.1 
Attachment resistance, RÂ (2) 1x 1075 
Bus resistance, Ri (Q cm7!) 4x 1075 

Cathode 
Conductivity pre-factor, A¢ (S Kem7!) 4.2 x 10° 
Conductivity activation energy, E¢ (J mol!) 9976.8 
Attachment resistance, RA (Q) 1x 10> 
Bus resistance, R, (2cm!) 8 x 1075 
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Fig. 3. Composite solution profiles for the base case. 


750° C and is heated by convection from the wall. However, as 
the flow proceeds along the channel, fuel is depleted and cur- 
rent density decreases. This, in turn, reduces the heat generation 
within the wall. Also, the tube is losing heat convectively to the 
exterior relatively cool air at 750°C. In the downstream sections 
of the tube, the flow carries energy from the upstream high- 
temperature regions downstream. In the downstream sections, 
the flow temperature slightly exceeds the tube-wall temperature. 
The velocity profile shows very little variation. This is because 
there is no molar change associated with oxidizing the syngas 
mixture. The shape of the velocity profile is caused primarily by 
the density variation associated with the temperature profile. 
Species profiles along the tube length show expected 
behavior. Hydrogen levels decrease owing primarily to electro- 
chemical oxidation. Carbon monoxide levels decrease owing 
primarily to water-gas-shift processes that react CO and steam 
to produce H? and CO2. The product species, H20 and CO2 
increase steadily along the length of the tube. The upper panels 
of Fig. 3 show gas-phase species profiles through the thickness 
of the porous anode at three locations along the tube. Again, 
expected behavior is observed. In the entry regions, strong gra- 
dients in hydrogen and steam are caused by charge-transfer 


reactions in the three-phase region near the dense electrolyte 
(top of the plots). Curvature in the profiles is caused by catalytic 
chemistry (primarily water-gas-shift) within the anode. In the 
downstream regions the gradients are diminished as the fuel is 
depleted and the current density is reduced. 

Fig. 4 shows the electric-potential profiles within the anode, 
anode bus, cathode, and cathode bus along the tube length. 
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Fig. 4. Local electric-potential profiles for the base case. 
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Fig. 5. Composite solution profiles for the reduced-attachment case. 


The terminal potential (at x = 0) is fixed at 0.65 V. Because 
of relatively low resistance in the axial buses, the electric- 
potential profiles along the busses are smooth. Also, because 
the attachment resistances are low in this example, there is very 
little electric potential difference across the attachment points. 
Although the axial busses have relatively lower resistance than 
the electrodes, they still have resistance. Thus, the cell potential 
(difference between cathode and anode potentials) must increase 
along the tube length, driving current toward the terminals at the 
left-hand side of Fig. 4. The current collection in this exam- 
ple is reasonably effective, with a net variation of cell potential 
between the axial buses of only about 0.05 V. 

It is evident from Fig. 4 that there can be significant electric- 
potential variation within the electrodes between the attachment 
points. Because of relatively high resistance within the elec- 
trodes, a potential difference is required to drive electric current 
axially along electrodes toward the attachment points. With the 
parameters in this example, local potential variations within the 
anode are about 0.015 V and within the cathode are about 0.03 V. 
Although the terminal potential is 0.65 V. the local cell potentials 
can be as high as 0.76 V in the entry regions of the tube (i.e., 
left-hand regions of Fig. 4). The locally higher cell potentials 
cause locally lower current densities. To maintain compara- 
ble electric-potential variations within both electrodes, many 
more attachment points are needed on the cathode, which has 
higher electrical-resistance. In the downstream regions of the 
tube, where the current density is diminished, the local electric- 
potential variations are also diminished. 

Figs. 5 and 6 show the model results for the reduced- 
attachment, where cathode attachment is reduced from 50 to 
25 points. This also causes the bus resistance between attach- 
ment points to double. Compared to the base case, reducing 
the cathode attachment causes much higher electric-potential 
variations between attachment points in the bus and within the 


cathode. Further, the reduced-attachment causes much greater 
variations in local cell potential and current density. Overall the 
effective operating potential of the reduced-attachment cell is 
substantially higher than the terminal potential. 

The local variations in current density are sufficient to cause 
noticeable variations in the MEA temperature (lower panel of 
Fig. 5). This is the result of locally varying heat generation 
within the MEA. As should be expected because of the con- 
vective nature of the flow, the temperature profiles within the 
gas flow are smooth. 

Because the local cell potentials across the MEA are higher 
for the reduced-attachment case, the local current densities are 
lower. This, in turn, causes lower heat release and thus lower 
cell and flow temperatures. Compared to the base case, the fuel- 
consumption rate is lower in the upstream regions and greater 
in the downstream regions. Despite the local high-frequency 
variations, the overall more-uniform fuel consumption rates and 
cell temperatures could be viewed as having some benefits. 

Fig. 7 shows the model results for the case without anode 
attachment. For anode-supported tubes with cathodes on the 
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Fig. 6. Local electric-potential profiles for the reduced-attachment case. 
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Fig. 7. Solution profiles without external anode current collection. 


outside, providing cathode current collection is generally easier 
than anode current collection where attachments must be made 
inside the tube. Thus, it is interesting to understand the impact of 
eliminating the anode current collection. This example assumes 
perfect cathode current collection (i.e., uniform cathode electric 
potential along length of the tube). With anode current collec- 
tion, all electric current must flow axially through the anode 
structure toward the terminal at the end. Therefore, the cell 
potential must increase along the tube length. For the parame- 
ters in this example, the cell voltage increases from the terminal 
potential of 0.65 V (at x = 0 cm.) to the open circuit potential of 
about 0.96 V within the first 10 cm of the tube. Once the open- 
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circuit potential is reached no further fuel can be consumed, and 
thus no more current can be produced. 

Eliminating anode current collection results in very poor per- 
formance for this example. The net power is low, and it is all 
produced within the first 10cm of the tube. Fuel utilization is 
low, with only a small fraction of the Hz consumed. Signifi- 
cant axial temperature gradients in the entry regions are also 
potentially problematic. 

Fig. 8 shows solution profiles for the case of uniform cell 
potential at Ece = 0.65 V. This is a limiting case, assuming 
perfect current collection on both the anode and cathode sides. 
Comparing Fig. 8 with Figs. 3 and 5, it is evident that the 
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Fig. 8. Solution profiles at uniform cell potential of 0.65 V. 
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uniform-potential case is more similar to the base case than the 
reduced-attachment case. 

It is interesting to note the tube-wall temperature variations 
in the vicinity of the tube inlet. The wall temperature at the inlet 
end is about 812° C for the base case, 790° C for the reduced- 
attachment case, 808 ° C for the without-anode-attachment case, 
and 840° C for the fixed-cell-potential case. These temperature 
variations, which are due to local heat release associated with 
chemical and electrochemical reactions within the anode and 
the ohmic losses across the MEA structure, are strongly coupled 
with the local current density. The local current density, in turn, 
is a strong function of the local temperature and the local cell 
potential. Therefore, it may be possible to design the current- 
collection system to control the cell-potential profile, seeking to 
achieve more uniform distributions of the current density and 
the temperature. 

Overall cell performance can be measured in terms of con- 
version efficiency, fuel utilization, and power density [29]. The 
cell efficiency is defined as 


22 We = JicFoen dA (36) 
Qin  MfinAhfin 


where We is the electrical work output and Qin is the heat that 
would be released upon full oxidation of the inlet fuel stream. 
The inlet fuel mass flow rate is min and Ahgin is the specific 
enthalpy associated with completely oxidizing the fuel stream. 
The electrical work is the product of the current density i, and 
operating voltage Ecen, integrated over the active membrane- 
electrode assembly (MEA) area. Fuel utilization Ucan be written 
as 


Mf out AN out 


U=1- (37) 


Mein Ahein 


where the “in” and “out” refer to the inlet and outlet of the fuel 
cell. The Ah refers to the specific enthalpy associated complete 
oxidation of any available fuels. This definition accounts for the 
energy content of any remaining fuels (or fuel byproducts) that 
leave in the fuel-cell exhaust. These definitions consider only 
performance within the SOFC, not overall system performance. 

Table 3 compares the cell performance from all cases. Except 
the case without anode current collection, the overall efficiency 
and fuel utilization for the other cases are comparable. The 
uniform-cell-potential case shows about 5% lower cell effi- 
ciency than the cases with discrete current collection. Since a 
large amount of fuel can not be consumed for the case without 
anode current collection, fuel utilization and cell efficiency is 
dramatically lower. 


Table 3 

Predicted overall performance for alternative current collection 

Case name Efficiency (%) Utilization (%) Power (W) 
Base case 52.5 95.3 32.4 
Reduced cathode attachment 52.9 90.4 32.7 
Without anode collection 14.6 24.1 9.0 
Uniform cell potential 48.6 97.2 30.0 


It is interesting to note that the efficiency and power den- 
sity are slightly higher in the reduced-attachment case than in 
the base case. Efficiency and power density depend on oper- 
ating potential, with efficiency usually maximum at operating 
potentials near 0.8 V [29]. The power density depends much 
more strongly on details of MEA performance, but usually 
peaks at lower operating potentials in the range of 0.6-0.7 V. 
The reduced-attachment case has a higher effective operat- 
ing potential, which tends to increase efficiency. However, the 
reduced-attachment case has lower fuel utilization. Thus, it 
should be possible to increase its performance by reducing 
fuel flow rate. The uniform-cell-potential case shows lower effi- 
ciency than either of the discrete current-collection cases. This 
is primarily due to the low operating potential of 0.65 V. A uni- 
form cell potential of around 0.8 V would lead to the highest 
efficiency. 


6. Summary and conclusions 


This paper develops and applies a computational model 
to assist current-collection design for tubular anode-supported 
solid-oxide fuel cells. The fuel cell model considers mass, 
momentum, and energy balances within the tube flow, reactive 
porous-media transport within the electrodes, electrochemical 
charge-transfer, and energy balances within the tube wall (i.e., 
the MEA). The current-collection electrical network considers 
multiple discrete attachment points to the cathode and anodes, 
as well as an axial bus that feeds current to the cell terminals. 

The primary purpose of this paper is to document the mod- 
eling approach, not specifically optimize current collection 
for a particular cell. Nevertheless, even with a few examples, 
some interesting observations can be made. Cell performance 
depends on the current-collection design. The resistances and 
the attachment points can significantly affect local temper- 
ature and species profiles along the tube length, as well 
as overall performance measures like efficiency, fuel utiliza- 
tion, and power density. Predictive models can be a valuable 
guide in the design of tubular SOFC cells and the required 
current-collection strategy. Certainly increasing the number of 
current-collection points leads to more uniform cell potential 
along the length of the tube. Because of higher electrical- 
resistance within typical cathodes, more attachments points are 
usually needed on the cathode than on the anode. However, 
even in cases with attachment separations leading to as much as 
150 mV variations in the local cell potential, overall performance 
measures such as efficiency and utilization are not strongly 
affected. 

The nominal cell used for illustration in this paper has a rel- 
atively small inner tube diameter (8 mm) and is relatively long 
(25 cm). This cell can generate net power around 30 W. With an 
operating potential of around 0.7 V, the axial buses must carry 
over 40 Amps. For the conditions used as examples in this paper, 
an anode bus made of silver needs to be over 4mm in diame- 
ter. Clearly a 4-mm-diameter wire within an 8-mm-diameter 
tube would significantly restrict flow area. The size and design 
of current collection, particularly inside anode tubes, can be a 
significant challenge in scaling to large tube dimensions. 
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